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■ Abstract 

O ! 

QQ , The yrast line of the N=Z nucleus °Cr is studied up to high spins by means of the 

0^ ' cranked Hartree-Fock-Bogoliubov method including the T=0 and T=l isospin pairing 

channels. A Skyrme force is used in the mean-field channel together with a zero-range 
density-dependent interaction in the pairing channels. The extensions of the method 
CJ \ needed to incorporate the neutron-proton pairing are summarized. The T=0 pairing 

^ ■ correlations are found to play a decisive role for deformation properties and excitation 

energies above IQfi which is the maximum spin that can be obtained in the f7/2 subshell. 

^ ' Pairing correlations play a crucial role in nuclear structure theory. They manifest themselves 
by the presence of a superfiuid condensate in the ground state of even-even nuclei, which is 
responsible for the odd-even staggering of many nuclear properties. In most nuclei, pairs are 
formed by identical particles (protons or neutrons) moving in time-reversed orbits[0]. However, 
one expects neutron-proton (np) pairing correlations to become important close to the N = Z 
line, since then protons and neutrons occupy the same shell model orbits[0]. 

Pairs formed by a neutron and a proton can be characterized by an isospin quantum number 
T equal either to 1 (as for pairing between identical nucleons) or 0. The properties of T=0 
correlations are by far not as well understood as the T=l correlations. In particular, not much 
is known about collective effects generated by the formation of T=0 pairs. 

Early studies of np correlations were mainly limited to light mass nuclei and schematic 
forces[^ ^, ^, 1^. The interest in np pairing correlations has been renewed^, ^ by the experi- 
mental possibilities opened by the event of powerful detector arrays. It is indeed now possible to 
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extend our knowledge on the structure oi N = Z nuclei to medium mass nuclei. In this paper, 
we will study the nucleus '^^Cr that has been investigated at Chalk River p and Legnaro |^ 
up to high angular momenta. 



The strongest evidence for np pairing comes up to now from masses. It has been shown |T0 



that the systematics of the nuclear masses for constant nucleon number A presents a cusp at 



N=Z and that this feature may be at least partly attributed to the np pairing[12, m. The 



T=0 pairs should also manifest themselves in the structure of excited states. The properties of 
rotational bands which have been an important tool to study the effects of particle-like pairing 
may also be appropriate for the study of np pairing. Cranked mean-field and Monte-Carlo 
shell-model calculations p, ^, |ll] have shown that the T=0 pairing field is more resistant 
against nuclear rotation than the T=l pairing. High spin states of N=Z nuclei may therefore 
be sensitive to the properties of T=0 pairing. 

In this letter, we present a study of the structure of high angular momentum states of ^^Cr. 
We have developed a cranked Hartree-Fock-Bogoliubov (CHFB) method with a Skyrme force in 
the mean-field channel and density- dependent zero-range pairing forces for the T=0 and T=l 
pairing channels. This is the first time that a Skyrme force is used in a calculation including 
the T=0 pairing channel. 

The yrast line of ^^Cr is known experimentally |^ up to an angular momentum / of 16h. 
It reveals a backbend in a plot of spin as a function of rotational frequency at a spin around 
12h. This backbend is related to a change from collective to non-collective states induced 
by the alignment of the f7/2 orbits|T^. The highest observed state has a spin of 16h which 
exhausts the f7/2 shell. A previous CHFB calculation |jl5| has shown that the ground state of 
^^Cr is prolate and that the shape evolves along the yrast line up to sphericity at a spin of 16h. 
"^^Cr has also been studied by the shell model |T^, [l^. Both shell model and CHF calculations 



reproduce successfully the backbending of the yrast line, although the CHF backbending curve 
is too low in angular velocity when compared to experiment. States with / >16h have not been 
determined in these works. The full fp shell has been included in the shell model space, but 
states with spins higher than 20h need an extension of the space to the sd-shell and the gg/2 
subshell. Such calculations are beyond present shell-model techniques. 

Since the mean field method describes rather successfully the low spin spectrum of ^*Cr 
without introducing T=0 correlations, we shall focus in this paper on the states beyond the 
'terminating' point at I=lQh. Our main aim is to determine to what extent the T=0 correlations 
will affect their structure. 

A new method of solution of the cranked Skyrme-HFB equations has been presented in 



recent publications [[17|, |T8[. It combines the imaginary-time evolution method to determine 
the basis which diagonalizes the mean-field hamiltonian and a diagonalization of the HFB 
hamiltonian matrix calculated in the HF basis to construct the canonical basis. Different 
interactions are used in the mean-field and pairing channels: a Skyrme force for the mean-field 



and either a seniority||T^ or a zero-range density-dependent interaction ||T^ for the pairing field. 
The pairing correlations were limited to like particle interactions. One of the important features 
of this method is a separation of the mean-field and pairing field equations, with however full 
self-consistency between them. 

To allow a simultaneous treatment of all forms of pairing, one has to generalize the Bogoli- 
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ubov transformation. As a consequence, the mean-field equations have to be modified: when 
pairing is hmited to hke particles, pairs are formed by states of opposite signatures. Since 
np pairs can occupy identical states, the signature is lost as a good quantum number by the 

Let and cf be respectively the neutron and proton 



IS 



general Bogoliubov transformation] 
annihilation operators in the states i diagonalizing the HF equations, with Ain neutron and 
Aip proton basis states. The Bogoliubov transformation from this single-particle basis to the 
quasiparticle basis in which the CHFB Hamiltonian is diagonal can be written as 



(c- \ 



V j 



(1) 



where a stands for an annihilation operator of a quasiparticle state, U ^ (^n^) are block ma- 
trices, and T denotes the transposition of a matrix. Note that the number of the quasiparticle 
operators {a} is A^n+A^p, and the size of U ^ and V ^ is (Afn+Alp) x M.^. JJ ^ and are 
(A4n+Aip) X Aip matrices. The quasiparticle states are now mixing protons and neutrons. 

It has been shown that when parity and signature are good quantum numbers and with 
an appropriate choice of phase, the HF wave functions can be represented in the x,y,z > 
octant of a box in r-space The complete description of the single particle wave functions 
requires then four real functions corresponding to the real and imaginary parts of the spin-up 
and spin-down components. The signature symmetry being lost, the space to be considered 
can only be reduced to a quarter of a box. 

m 



The method introduced in ref 



can be rather easily extended to this new case. After 
diagonalization of the Bogoliubov equations written in the HF basis, the density matrix can be 
constructed and diagonalized. The wave-functions in the canonical basis have now 8 components 
since proton and neutron states are mixed by the pairing correlations. The density is given by: 



p(r) = J2v^\ipi 



(2) 



where <fi{r,a,T) is the complex (cr, r) component of the i^^ canonical basis wave-function; a 
and T denote the spin and isospin respectively and vf is an eigenvalue of the density matrix. 

Parlihska et al. [^] have derived the new terms that arise in a Skyrme force when the 
densities mix protons and neutrons. For simplicity, and since we use different interactions in 
the mean-field and pairing channels, we have not introduced these terms. Densities for protons 
and neutrons are defined by limiting the summations in equ. (2) to r equal either to 1/2 or 
— 1/2. The occupations of neutron and proton states is then given by / (PrJ2i,cr'^i\'^i{'^^ 

In the following, we have used the parameter set SHI for the Skyrme interaction |p2l and a 



zero-range pairing force as in ref. [[18 



VTiri, ra) = 6*^(1 - p(ri)/pc) 5(ri - ra) 



(3) 



where Gt is the strength of the interaction and pc a constant taken equal to 0.16 fm ^, a 
value close to the nuclear saturation density. In this way, Vr(ri,r2) is peaked at the surface 



3 



of the nucleus. Assuming isospin invariance, the same strength has been used for all the T=l 
components of the pairing force. We have fixed Gt=i to 1000 MeVfm^, which is the value 
used in a previous HFB calculation of the ground state properties of Mg isotopes far from 
stability |2^. As usual, a smooth cut-off of the pairing force has been introduced 5 MeV above 



the energy of the Fermi level. For the strength of the T=0 correlations, we have chosen to 
scale the T=l strength by a factor equal to 1.3. This choice originates from proton-neutron 
scattering data that show clearly the spin singlet (>S'=0) interaction to be weaker than the 



triplet (5=1) one, resulting in a ratio of ~ 1.3 for Gt=o/Gt=i [0 



Using the antisymmetry of the pairing tensor k and with a calculation similar to that of 



ref|18|, one can write the pairing field A in a way suitable for computation and decompose it 



into its T=l and components: 



era' 



11 1 1 



kl (Tcr 

1 1 „ „. , , 1 . , , 1 



At=i,tM = -i;T=i(r)^^(-r-r'|lT^)0fc(r,-,r)0i(r,--,r')K(fer)(Zr') , (6) 



kl 



TT' 



VT=i,o{r) = Gt=i,o (1 - p(r)/Pc) , (7) 

where 0j(r,(T, r) is the i^^ single-particle wave function in the HF basis. 

The spin singlet S=0 with Tz=±l term in equ. (6) corresponds to the usual pairing inter- 
action between signature-reversed states. The np pairing interaction appears in both the spin 
singlet (T=l) and the spin triplet {T=0) modes. Note the symmetry between the spin anti- 
symmetric (isospin symmetric) state (S'=0, T=l) and spin symmetric (isospin antisymmetric) 
state (5=1, r=0). 

With the symmetries imposed in previous works |^ on the single-particle wave-functions, 
a real pairing tensor k does not allow the simultaneous presence of all the pairing modes. A 
complex K is then generated by a complex Bogoliubov transformation. The situation is different 
when, as here, parity is the only symmetry and the sp wave-functions are complex. When the 
system is time-reversal invariant, one can easily show that T=0 and T=l pairing may be present 
at the same time. Using the relation between time reversed states given in ref one obtains 
a purely imaginary AT=o,5,=o(r): 

ImAT=o,s.=o(r) 

= -V2vT=oir)Y,i-)^''^Y. \ Im0j(r, -(T, -r)/t(i^)(j_^) 

+ X^Re0,(r, or, r) Im0j(r, -a, -r)K(.^)(j_^) \ , (8) 

ij ) 
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Figure 1: Pairing energies corresponding to the four pairing modes, along the yrast hne. The energy 
of the T=l, Tz = (np) component is represented by diamonds. The inset shows the expectation 
values of the total intrinsic spin along the rotational axis at high angular momenta. Squares and 
triangles correspond to positive and negative parity, respectively. 

where j denotes a time- reversed state of j. In the same way, one can show that Ar=i 5^=0(1") 
is reaL Since the components of the pairing potential for all values of Sz and may be 
present at the same time with a real pairing tensor, we have considered only real Bogoliubov 
transformations. In most of our results however, our solutions are dominated by either the 
T=0 or the T=l pairing field. This difficulty of obtaining mixed solutions has been discussed 
in ref. and is due to the lack of particle number and isospin conservations. 

With a real pairing tensor, one can show that both the pairing and mean-field matrix 
elements are real, even when time reversal invariance is not assumed and signature is not a 
good quantum number. A detailed analysis of our method will be presented in a regular paper. 

At low spins, we have found two solutions dominated by either the T=l and T=0 pairing 
modes. Up to a spin of 16 h, the two bands differ by their pairing properties but also by their 
deformation and excitation energies. For spins between 4 and 8 h, the T=0 solution displays 
both pairing modes, illustrating the ability of our method to describe them simultaneously. We 
defer a detailed discussion of its properties to a forthcoming paper and limit our discussion 
of the low spin spectrum to the T=l solution. The evolution of the four components of the 
pairing energy as a function of the angular momentum is plotted on Fig. 1. Pairing correla- 
tions disappear rapidly and for angular momenta between 6 and 16 h, a pure CHF solution 
is obtained. At 16^, one reaches the total alignment of protons and neutrons in the f7/2-shell. 
The higher spin states involve particle-hole excitations from the sd- and f7/2-shells into the gg/2- 
and f5/2-shells. Above 1=16^, the T=0 pairing channel becomes active and persists up to the 
highest spins that we have studied. These T=0 correlations reach a maximum at / = 24h. 

The calculated spins as a function of frequency are compared to the experimental data 
on Fig. 2a. Up to / = IQh, our results agree with previous CHFB calculations using the 
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Figure 2: a)Variation of the angular momentum as a function of the rotational frequency along the 
yrast line. Our results are represented by □ for Gt=o = 1.1Gt=i and A for Gt=o = 1-3Gt=i in the 
upper part. The experimental data (•) are taken from ref.|^, b)The contributions to the angular 
momentum coming from positive (squares) and negative parity states (triangles). 




Figure 3: Variation of the quadrupole deformation (parameterized by /3 (circles) and 7 (squares)) as 
a function of the rotational frequency along the yrast line. 
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Gogny force. Like in those calculations, the spectrum is slightly too compressed with respect to 
the experimental data but the general agreement is satisfactory. Since pairing correlations are 
weak at low spin in '^^Cr, one can expect particle number projection techniques to decrease the 



moment of inertia at the bottom of the band and to improve the agreement with the dataflS 

Above 16 h, the structure of the yrast line changes. The 7 transition energies become almost 
constant, around 5MeV, with two slight shifts at 24 and 32h. As illustrated on Fig. 2a, the 
results are not affected by decreasing the pairing strength from Gt=q = I.SGt=i to 1.1Gt=i- 
Note that solutions in the spin region / = 18 — 30h are obtained only in the presence of T=0 
pairing correlations. If this mode is switched off, we have not been able to find neither CHFB 
nor pure CHF solutions; states can only be constructed as non collective multi-qp (or sp) states. 

To generate spins above 16^, one can either promote particles from the subshell into 
the fp shell above the N(Z)=28 gap or via particle-hole excitations from 0^3/2 into gg/2- For 
both protons and neutrons, excitations into the yield a maximum gain in spin of 2h while 
the (^3/2)"^ —99/2 ph excitations generate increases in spin up to I=6h. At /=0, the energy of 
this 2huj ph-excitation is of the order of 14 MeV but it is drastically reduced by rotation. For a 
cranking frequency around 2.5huj MeV, the 1^3/2 and gQ/2 levels are close to the Fermi surface. 

The contributions of the different parities to the total spin are depicted in Fig. 2b. The 
main contributions above 16h are coming from positive parity states. In our calculations, the 
occupation of the gg/2 {(^3/2) grows (decreases) monotonically from 0(1) to 1(0.4) in the spin 
range of 16h to 32h. It reaches 0.6(0.73) at / = 24h where the T=0 pairing energy has its 
maximum. However, np-pairs of negative parity contribute also to the T=0 correlations. At 
32h, although pairing correlations are still sizable, the angular momentum is mainly generated 
by a few single particle levels: (g9/2 dj^/2)6+ ® (f7/2)^6+ ® ih/2 ^5/2)4+ for both neutrons and 
protons. 

The evolution of the quadrupole deformation calculated along the yrast line is shown on 
Fig. 3. At low spin, "^^Cr is prolate; its deformation decreases above I = Ah and it becomes 
spherical at / = 16^. Except for the small increase of deformation at low spin, our result is 



very similar to the one obtained with the Gogny force[15|. Above / = 16/;., the deformation 



is increasing smoothly with angular momentum. The nucleus is calculated to have an oblate 
shape rotating around the symmetry axis. 

The deformation energy surface at / = 24/i is shown on Fig. 4. It presents only an oblate 
minimum rather soft against deformation. Q This may indicate that one should include dy- 
namical fluctuations along the 7 degree of freedom. 

It is well established that deformation may be induced by the long range quadrupole- 
quadrupole force due to particle-hole excitations like ^3/2 — 5'9/2- The T=0 correlations induce a 
different and new type of collective excitations. The quadrupole deformation is now enhanced 
by the short-range T=0 np-force. The magnitude of the T=0 pairing correlations reflect the 
pair scattering from a np ^3/2 to a np (79/2 pair. Since these pairs carry alignment, the smooth 
increase of their occupation result in an increase of the angular momentum which is of course 
correlated with the increase in deformation. Spin is also generated by the scattering of a np 
pair from the f-j/2 to the f^/2 orbital. Both processes are only possible because the T=0 mode is 

"'^We also performed cranked Strutinsky type calculations, that show the existence of oblate super-deformed 
and prolate hyper- deformed solutions at « lOMeV and ~ 13MeV above the yrast /=16?i state. 
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Figure 4: Potential energy surface calculated at / = 24h. The minimum of the energy is obtained at 
/3 = 0.11, 7 = 60°. 

present. Thus, T=0 pairing correlations provide a new mechanism to generate angular momen- 
tum and quadrupole deformation beyond 'terminating states'. Note also that the T=0 pairing 
correlations become active in a spin regime, where one does not expect pairing to play any role. 

The expectation value of the total intrinsic spin < 5*^ > gives additional informations on the 
structure of the np pairs. Since the ^3/2 pair has its intrinsic spin < Sz > coupled anti-parallel 
to the orbital angular momentum, the scattering into gQ/2 results in an increase of < 5"^ >. The 
opposite holds for the scattering from f'j/2 into The increase (decrease) in < Sz > for the 
positive (negative) parity is shown in the inset of Fig. 1. Interestingly, in algebraic models, the 



main building blocks for neutron-proton pairing]^ are constructed by the coupling of a pair 
of spin orbit partners, like /5/2 and /7/2 to L = and S = J = 1. In general, this couphng 
scheme is quenched above A=40 because the spin-orbit splitting becomes too large. We suggest 
that the reduction in the Wigner energy when going from the sd-shell to the fp-shell noticed 



by Satula et al. |13[, may be traced to this effect. We find here that such spin-orbit pairs may 
still play a role at high angular momenta. 

A pure HF state can be constructed at I=32h. Its deformation properties are compared to 
the CHFB solution in Table 1. Although pairing correlations are weak at I=32h, the spreading 
in occupation probabilities that they generate results in a quadrupole moment much larger 
than that of the CHF state. On the other hand, the g-factor of the two solutions is quite close, 
around 0.52. In fact, the g-factor is nearly constant all along the yrast line and close to the 
value obtained with a Gogny force ||T3|. 

In summary, we have studied the properties of the high angular momentum states of ^^Cr, 
beyond the 'band termination' by means of the CHFB theory with both the T=0 and T=l 
pairing channels taken into account. The only symmetry imposed on the wave-functions permit 
to limit oneselves to real Bogoliubov transformations and hence to have a real pairing tensor 
K. We have shown that, contrary to previous studies where complex-HFB transformations had 
to be introduced, the two types of pairing interactions can be simultaneously present. 
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CHF 


CHFB with r=0 pairing 


Qo 


total [fm^] 


123 


205 


Qc 


charge [efm^] 


63 


104 






0.16 


0.26 


7 


[deg] 


60 


59 



Table 1: Total quadrupole moment (Qo), charge quadrupole moment (Qc) and deformation P 
and 7 at I—32h. 7=60° means oblate noncoUective state. 

The yrast hne above 16h has properties intimately linked to the presence of T^O pairing 
correlations: smooth increase of the quadrupole moment as a function of spin and nearly 
constant transition energies in the range 4.8 to 5.6 MeV. The T=0 correlations maximize 
exactly at half the maximum spin value of the involved configurations and induce a new type 
of collectivity. A spectacular effect of the np T—0 correlations is the presence of a much larger 
deformation at / = 32h than in the pure CHF state. All these properties should be fingerprints 
of the presence of T=0 pairing at high spins. 
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